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ABSTRACT 

Supermassive black holes ejected from galaxy nuclei by gravitational wave recoil will carry a retinue 
of bound stars, even in the absence of an accretion disk. We discuss the observable signatures related 
to these stars, with an emphasis on electromagnetic flares from stars which are tidally disrupted by 
the black hole. We calculate disruption rates for the bound, and the unbound, stars. The rates are 
smaller than, but comparable to, rates for non-recoiling black holes. A key observational consequence 
is the existence of powerful off-nuclear and intergalactic X-ray flares. We also discuss other observable 
signatures associated with the bound stars, including episodic X-ray emission from accretion due to 
stellar mass loss, intergalactic supernovae, and feedback trails. 
Subject headings: galaxies: active - galaxies: evolution - galaxies: X-rays 



1. INTRODUCTION 

Gravitational waves emitted during gravitational col- 
lapse carry away linear momentum, causing the center 
of mass of the collapsing object to recoil (Peres 1962, 
Bekenstein 1973). In the case of binary black holes, coa- 
lescence can be accompanied by kicks as large as ~ 200 
km s^^ if the holes are nonspinning prior to coalescence, 
increasing to ^ 4000 kms~^ in the case of maximally 
spinning, equal-mass holes with the optimum orienta- 
tions (e.g., Campanelli ct al. 2007; Gonzalez et al. 2007; 
Herrmann et al. 2007; Koppitz et al. 2007; Baker et al. 
2008; Schnittman et al. 2008). 

Recently, the first compelling candidate for a recoil- 
ing SMBH was discovered (Komossa et al. 2008). The 
quasar SDSSJ092712. 65+294344.0 at z = 0.7 exhibits 
two separate sets of unusual emission lines: a set of broad 
emission lines, presumably associated with gas bound to 
a SMBH, and a second set of atypically narrow emission 
lines. The broad lines are shifted by ^ 2650 km s^^ rela- 
tive to the narrow lines. This is one of the key predicted 
signatures of a recoiling SMBH, which would carry with 
it the gas from the broad line region while leaving the 
bulk of the narrow line gas behind (Merritt et al. 2006; 
Bonning et al. 2007). 

Gas near to the SMBH is expected to respond to the 
kicks in a number of other, potentially observable ways. 
Accretion activity may be temporarily interrupted dur- 
ing the final phases of binary coalescence but would re- 
start thereafter with a certain time delay (Liu et al. 
2003; Milosavljevic & Phinney 2005). The recoiling 
SMBH would then appear, temporarily, as a "quasar" 
that is spatially offset from the core of its host galaxy 
(Madau & Quataert 2004, Loeb 2007). UV (Lippai et al. 
2008), soft X-ray (Shields & Bonning 2008) and IR flar- 
ing (Schnittman & Krolik 2008) could result from shocks 
in the accretion disk surrounding the SMBH just after 
recoil, or when the inner disk reforms. 

All of these detection methods require the presence of 



an accretion disk around the coalescing SMBHs, and the 
duration of the predicted signal is limited to the time it 
would take the bound gas to be accreted, or much less. 
But such gas may be lacking in mergers of early type 
galaxies, and even in gas-rich mergers, inspiral of the bi- 
nary SMBH may stall long enough that the nearby gas 
is mostly depleted. How can we observe such recoils ? 
Even in the absence of a gaseous accretion disk, a recoil- 
ing SMBH will always be accompanied by a retinue of 
tightly bound stars (Gualandris & Merritt 2008), and a 
fraction of these stars will eventually undergo tidal dis- 
ruption. Stellar tidal disruptions (Luminet 1989) have 
been observed in the form of giant-amplitude, luminous 
X-ray flares (Komossa & Bade 1999; Halpern et al. 2004; 
Komossa et al. 2004). In the case of a recoiling SMBH, 
these flares will appear off- nuclear. In this Letter we 
present tidal disruption of stars as one of the most uni- 
versal observational signatures of recoiling black holes, 
and discuss other signatures related to bound stars. 

2. DISRUPTION RATE OF BOUND STARS 

A recoiling SMBH carries with it a cloud of stars on 
bound orbits. We first compute the mass of this cloud 
(expressed as a fraction /f, of the SMBH's mass), then we 
compute the rate at which the stars would be scattered 
into the SMBH's tidal disruption sphere. 

Figure la shows the steady-state distribution of bound 
stars following an instantaneous kick of magnitude Vk, 
assuming an initial, power-law density profile p = 
Po{r/ro)^^ around the SMBH. ^ Stars initially at dis- 
tances r ^ rk = GMbk/V^ from the SMBH will be un- 
bound after the kick. For T4 ^ 0.4T4sc i-e. large enough 
to remove the SMBH from the galaxy core (Gualandris 
& Merritt 2008; Ksc is the escape velocity), ry^ is small 
compared to the SMBH influence radius and stars that 

^ The bulk of the recoil velocity is imparted to the co- 
alesced SMBH in a, time ~ GMbh/c^ (e.g. Fig. 1 of 
ICampanelli et al. 20071 ) and so the kick is essentially instantaneous 
as seen by the stars. 



Komossa & Merritt 



follow the SMBH will be moving on essentially Keplerian 
orbits both before and after the kick. The final distribu- 
tion can be computed uniquely from the initial distribu- 
tion by randomizing the orbital phases of the stars that 
remain bound. Beyond ^ r^ the cloud consists of an 
elongated, p ~ r~^ envelope containing stars that were 
pushed into nearly unbound orbits by the kick (Fig. lb) . 
The bound mass must scale as Mb ~ p(''k)''k ~ 
PQrl{GM^Y{/rQV^Y^'^ . Choosing for tq the (pre-kick) 
influence radius r, of the SMBH, defined as the radius 
containing a mass in stars equal to twice Mbh, this be- 
comes 
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BH 
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we find F(7) « 11.67"i'^'5, 0.5 ^ 7 ^ 2.5. An alternative 
form is /,,_3 = G(7)A<7 r;^nfl'iK'ir'^(^k/K,sc)^(''-3) 
where /b^^g = ^/10-^ M,j = Mbh/IO^M©, n, 



'infl,l 



r./lO pc, Ksc,3 = 14sc/10^ km s'^, and 6(7) = 
(0.048, 0.22, 1.6, 15) for 7 = (0.5, 1, 1.5, 2). For kick ve- 
locities in the range of interest, equation ([1]) predicts 
bound masses of order 1% of the SMBH mass, falling 

as V^ '^ . The stars around a recoiling SMBH would 
resemble a globular cluster in total luminosity, but with 
a much greater velocity dispersion due to the large bind- 
ing mass Mbh; observationally they might look like ul- 
tracompact dwarf galaxies (e.g., Drinkwater et al. 2003). 
The appropriate value for 7 is the density slope just 
prior to the kick, and after the binary SMBH has com- 
pleted its inspiral; r, is likewise defined at this time. 
Slow inspiral of the SMBHs in a large, low-density galaxy 
produces a flat core, 7 w 0.5 (e.g. Merritt 2006); how- 
ever it is not clear whether binary SMBHs in such nuclei 
can overcome the "flnal parsec problem" (Milosavljevic 

6 Merritt 2003) and coalesce. In lower-mass galaxies, 
Mgai ^ 10^°A^o, binary evolution can continue to coa- 
lescence in < 10 Gyr via coUisional loss-cone repopula- 
tion; in this case the pre-kick density profile would have 

7 ~ 1.75 (Merritt et al. 2007). Rapid, gas-driven inspi- 
ral would tend to preserve the initial density profile and 
might even steepen it via star formation. Henceforth we 
take 7 = 1 as a fiducial value. 

In standard loss-cone theory, stars are scattered into 
the tidal disruption sphere r < rt at a rate ^ 
N(r)/ [tr{r)l'a{r/rt)] where N{r) is the number of stars 
within r and tr{r) is the local (non-resonant) relaxation 
time.^ The total disruption rate from the cloud of bound 
stars would be 



Nnr ~ Cnr{i) 



InA /Vk 
ln(rk/rt) \r]^ 



fl 



(2) 



with InA « ln(MBH/™^) the Coulomb logarithm; the 
ratio of logarithmic terms is of order unity. Equa- 
tion (m assumes that gravitational encounters are un- 
correlated; however, near the SMBH, orbits are slowly 
precessing Keplerian ellipses and "encounters" are highly 
correlated (Ranch & Tremaine 1996), shortening the ef- 
fective, angular momentum relaxation time by a factor 
^ m^N [r) / M-Q-R (Hopman & Alexander 2006). The con- 

^ The logarithmic factor is associated with diffusive loss-cone 
repopulation which is appropriate in the "empty loss cone" regime 
near a SMBH (e.g. Lightman & Shapiro 1977). 



tribution of this "resonant relaxation" to tidal disrup- 
tion rates of SMBHs embedded in galaxies is small since 
most of the disrupted stars come from orbits with r Ki r,. 
However in our case, stars beyond ^ rk <IC r, were re- 
moved by the kick and loss-cone repopulation is domi- 
nated by resonant scattering, yielding 



Nrr ~ Crr{j) 



InA 



ln(rk/rt) V rk 



^1/. 



(3) 



i.e. Nun > N^n for A < 1. 

The coefficient Crr in equation ^ is poorly deter- 
mined (Ranch & Ingalls 1998; Hopman & Alexander 
2006). Because the number Alh/m^ of stars remaining 
bound to a recoiling SMBH is relatively small, tidal dis- 
ruption rates in this regime can be computed directly 
via sufficiently accurate TV-body integrations. Figure Ic 
shows the results of a series of such experiments starting 
from initial conditions generated from the steady-state 
(nonspherical and anisotropic) distribution of Figure la, 
realized with TV = 1.5 x 10^ particles and various values 
of m ^. We used t he hybrid TV-body code (/)GRAPEch 
(Harfst et al. 20081) running on the RIT GRAPE cluster 
(JHarfst et al. 2007 ). Stars were initially removed from 



the model if their periastra fell below 10^ rk, the as- 
sumed radius of the disruption sphere; the model was 
then integrated forward and stars that approached the 
SMBH particle at distances < rt were recorded and re- 
moved from the integration. The TV-body integrations 
confirmed the TV oc /{, dependence predicted by equa- 
tion dSI). ^ In the adopted units (G = Mbh = Vk = 1), 
the disruption rate is ~ 0.15/;,; assuming that A w 
rk/(2GTOj,/y^) « AIbh/itt-* in the simulations then im- 
plies CflR(7 = l)« 0.14. 

Using the derived value of Cru, equation ([3]) predicts 
for the disruption rate of bound stars (7 = 1): 



TVb 



6.5 X 10-6yr-iAf.-;ivi3^/,^_3 



(4) 



here and below. In A/ln(rk/rt) has been set to 2. Com- 
bining equation ([4]) with equation ([T]) for the bound mass 
gives 

TVb « 1.5 X lO-'yr-'M,jrr^l,V,-^. (5) 

Figure 2 (solid lines) plots TVb for SMBHs ejected 
from the centers of two representative galaxies, with 
masses 4.5 x IO^^Mq and 1.5 x IO^Mq; the galaxies were 
modelled as Prugniel-Simien (1998) (i.e. de-projected 
Sersic) spheroids with Sersic indices (4,2.5). How do 
the rates we estimate compare with flare rates of non- 
recoiling SMBHs? The steady-state stellar disruption 
rates of SMBHs embedded in galactic nuclei are pre- 
dicted to be in the range 10~^yr~^ ^ TV ^ 10~^yr~-'^ for 
Mbh ~ lO^M© (e.g. Wang & Merritt 2004, Fig. 5b). 
For Vk ^ 10"^ km s~^, rates derived here are roughly an 
order of magnitude lower. Were resonant relaxation not 
effective, this ratio would be ^ 10"^ rather than ^ 10. 

Ignoring changes in the shape of the density profile 
of the bound population with time, the disruption rate 
is predicted to drop off as TV = TV(0) exp(— f/r), r w 

^ We note that the various other conditions required for resonant 
relaxation to be present in these simulations were satisfied, i.e. the 
integration times were long compared with orbital precession times 
and most of the stars were in the diffusive, as opposed to pinhole, 
loss-cone regime. 
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3.6GM|jj/l^^m+. This decay is included in the curves 
of Figure 2. In fact, after ~ 1 (non-resonant) relaxation 
time, the density would evolve to the Bahcall-Wolf p oc 
r ''^ slope and maintain this profile as its amplitude 
decayed; this complication was ignored in the curves of 
Fig. 2. 

We note that any post-merger galaxy would likely 
be highly inhomogeneous. Tidal disruption rates could 
be temporarily increased during close encounters of the 
SMBH to a massive perturber, e.g. a giant molecular 
cloud, in much the same way that "comet showers" are 
triggered b y near passage of the solar system to a star 
(|Hills 1981h . The same is true for the rate of disrup- 
tion of unbound stars (§3) if the SMBH passes through 
a dense clump. 

3. CONTRIBUTION OF UNBOUND STARS 

Even a naked SMBH encounters stars as it traverses a 
galaxy (Kapoor 1976). Unbound stars are deflected into 
the tidal disruption sphere at an instantaneous rate 

^l.7xlQ-'>yY-^~p{r/R,)p,^iMt'jV-^R,^ii (6) 

where p is the local (stellar) density, p = p/p{Re), Re is 
the galaxy effective (projected half-light) radius, pe,i = 
/9(i?e)/lM0pc-3, and R^^ii = RJlO^^ cm; V, = K(r) 
is the instantaneous SMBH velocity. Figure 2 (dashed 

lines) shows iVunb('') in two representative galaxy mod- 
els. (The curves in Fig. 2 include the correction factor of 
Danby & Camm (1957) that accounts for the decrease in 
the capture rate when V, « a, with a the stellar velocity 
dispersion.) In the larger of the two galaxy models con- 
sidered, of order 10 flaring events are predicted from both 
bound and unbound stars for Vk = 10^ km s^^during the 
time required for the SMBH to travel from the core to 
the half-mass radius; in the smaller galaxy, N^nh ^ -^b- 

4. OBSERVABILITY 

We have shown that rates of stellar disruption by re- 
coiling SMBHs are interestingly high; only moderately 
lower than those of SMBHs in the cores of galaxies. For 
a typical galaxy, we predict r^ 20 flares as the SMBH 
travels through the central parts of the galaxy (50% of 
these from bound stars, 50% from unbound stars), and 
~ 10^ - 10^ more after the SMBH has left the galaxy. 
Now, we discuss observational consequences. 

4.1. Powerful off-nuclear X-ray flares and feedback trails 

Stellar tidal disruptions appear as luminous X-ray 
flares which reach quasar-like luminosities and then de- 
cline on month-year time scales (e.g. Komossa & Bade 
1999). Tidal flares from recoihng SMBHs would look 
similar but would occur off-nucleus, and would be eas- 
ily identified since no other known mechanism produces 
quasar-like luminosities far from the nucleus. 

Future X-ray all-sky surveys like eROSITA (Predehl et 
al. 2006) and EXIST (Grindlay et al. 2005) will be sensi- 
tive to these types of events. While sky surveys are most 
suited to identifying the candidates, follow-up observa- 
tions with high spatial resolution (~ 0.5", Chandra) will 
then allow the measurement of spatial offsets from the 
galaxy core. If the flare is associated with jet emission, 



radio observations would greatly improve the positioning 
accuracy. With future spectroscopic instrumentation, we 
will be able to determine the line-of-sight recoil velocity 
directly, if emission lines form in the temporary accretion 
disk of the stellar debris. Such instrumentation will be 
available aboard XEUS, its NFI reaching 3 eV resolution 
at 6 keV (Turner et al. 2007). While amplitudes of vari- 
ability are highest in the soft X-rays, optical wide-field 
surveys like Pan-STARRS (Kaiser et al. 2002) will also 
be sensitive to flare events. 

Before the eROSITA mission, to be launched in 2011, 
will start its all-sky survey, during the first half year of its 
operation a deep survey of the nearest clusters of galax- 
ies will be performed (G. Hasinger, priv. com.). This 
will provide an excellent opportunity to detect the near- 
est flare events. Disruption rates of unbound stars are 
highest in the galaxy core (Fig. 2); in order to iden- 
tify these as off-nuclear events, high spatial resolution is 
required (at the distance of the Virgo cluster, 1" corre- 
sponds to ^100 pc). The long time scales for damping 
of orbital oscillations (Gualandris & Merritt 2008) when 
Vk < Vi^sc will increase the chances of detecting an offset. 
In addition, the abundance of "fossil" intracluster recoil- 
ing SMBHs from past mergers (Volonteri 2007) would be 
especially high in clusters. 

While the brightest phases of a tidal disruption flare 
might only last a few years or decades, interaction with 
the environment would leave a feedback trail along the 
path of the SMBH. If radio jets are temporarily formed 
during the accretion phase (Wong et al. 2007), they 
would create local cavities in the ISM, similar to but 
smaller than the X-ray cavities that have been observed 
in nearby ellipticals (e.g., Biller et al. 2004) and that 
have been suggested to be linked to tidal disruption 
(Wang & Hu 2005). The bright flare wiU also excite 
emission lines in any surrounding gas. Even though these 
lines would be faint given the short time span of the en- 
ergy input, low density gas would retain memory of the 
flare for a long time (e.g., for a particle density n = 10 
cm~^, the hydrogen recombination time scale is on the 
order of 10000 yrs). This way, quasar-like emission-line 
ratios could be produced far away from the nucleus of a 
galaxy. 

In addition to stellar disruption and related observa- 
tional signatures, mass loss from evolving bound stars 
(e.g., Kudritzki & Puis 2000) will provide episodic low- 
level fuelling of the SMBH, causing repeated episodes of 
X-ray emission and variability. 

4.2. Intergalactic flares and other signatures 

Tidal disruptions continue long after the recoiling 
SMBH (with Vk > V^sc) has left its host galaxy and wan- 
ders in intergalactic space (Fig. 2). Such events will 
appear as luminous flares without host galaxies. Fur- 
thermore, the bound stars will undergo stellar evolution 
and will therefore ultimately produce intergalactic plane- 
tary nebulae, and intergalactic supernovae (SNe) of type 
I '^. While the most massive recofling SMBHs {M > 10^ 

■* We note that high-mass stars that produce type II SNe evolve 
quickly; within t < 100 Myr, a SMBH moving at 1000 km s~i 
has reached 100 kpc. Therefore, most of the type II SNe would be 
produced when the SMBH has not yet escaped into intergalactic 
space. These events might be identified from their large velocities 
with respect to their host galaxy. 
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Mq) would generally not disrupt solar- type stars, they 
would still partially disrupt giant stars, and they would 
also show repeated phases of activity from fuelling due 
to stellar mass loss. As such they might hide among the 
"blank field sources" with bright X-ray emission but no 
optical counterpart (e.g. Cagnoni et al. 2002) if the ab- 
sence of optical hosts persists in deeper optical imaging. 
After several Gyrs, stars will start turning into white 
dwarfs, and SMBHs which left their galaxies long ago will 
have a cloud of white dwarfs, neutron stars and stellar 
mass BHs bound to them. While the white dwarf tidal 
radius is inside the Schwarzschild radius except for low 
( ^ lO^-^Afo) BH masses. Dearborn et al. (2005) have 
shown that relativistic compression causes white dwarfs 
to ignite and explode as SN at distances up to lOOi?^ 
from a SMBH. Part of the disrupted white dwarf will 



escape while the rest will eventually be accreted; both 
the SN signal and the accretion phase would allow us to 
detect such fossil intergalactic SMBHs. 

In summary, we have discussed observational signa- 
tures of recoiling black holes related to the bound (and 
unbound) stellar population. Finally, we emphasize that 
all these signals would generically be associated with 
recoiling SMBHs, whether or not an accretion disk is 
present initially, and that they would continue episodi- 
cally for a time of ^ 10 Gyr. 
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Fig. 1. — (a) Contours of the projected density of a cloud of bound stars around a kicked BH. x and y are spatial scales and the unit 
of length is GMbh/V^. The initial density was p oc r~^ and the kick was in the x-direction. Contours are separated by 0.17 in the log. 
(b) Spherically-symmetrized density profiles of stars around kicked BHs; the initial profiles — p oc r"'', 7 = (1, 1.5, 2) - are indicated by the 
dotted lines. Density normalization is arbitrary, (c) Dependence of disruption rate on bound mass in a set of Af-body integrations with 
N = 1.5 X 10* and various values of the mass Mj in stars bound to the BH. dN/dt is in A^-body units, as defined in the text. The dashed 
line has unit slope, the expected dependence if feeding is dominated by resonant relaxation. 
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Fig. 2. — Tidal disruption rate as a function of distance travelled for SMBHs ejected from the centers of two galaxy models. Model 
1 (black lines): Mg^i = 4.5 X IO^'^MqjMbh = 3 x 10^MQ,Sersic index n = 4.0, i?e = lkpc,r, = lOpc. Model 2 (red lines): Mg^i = 
1.5 X 1O9M0,Mbh = 1 X 10'5MQ,n = 2.5, iJe = O.Skpc, r. = 3pc. Line width indicates kick velocity: Vi = 500km s-i(thinnest), 
1000,2000,3000,4000 km s~^ (thickest). Solid lines show N for bound stars, dashed lines for unbound stars. Lines terminate at the right 
where the SMBH reaches apocenter; distances r <r, are indicated via dotted lines. The falloff in the bound rates at large radii is due to 
depiction of the bound population. The tick marks indicate time and are at log (t/yr) = 5,6,7,8,9. 



